Parallel Transport of Electrons in Graphene Parallels Gravity 



Andrej Mesaros^, Darius Sadri^^^, and Jan Zaanen^ 
^ Instituut-Lorentz for Theoretical Physics, Universiteit Leiden, 
P. O. Box 9506, 2300 R A Leiden, The Netherlands and 
'^Institute of Theoretical Physics, Ecole Polytechnique Federale de Lausanne (EPFL), CH-1015 Lausanne, Switzerland 

Geometrically a crystal containing dislocations and disclinations can be envisaged as a 'fixed 
frame' Cartan-Einstein space-time carrying torsion and curvature, respectively. We demonstrate 
that electrons in defected graphene are transported in the same way as fundamental Dirac fermions 
in a non-trivial 2-1-1 dimensional space-time, with the proviso that the graphene electrons remember 
the lattice constant through the valley quantum numbers. The extra 'valley holonomy' corresponds 
to modified Euclidean symmetry generators. 
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Introduction -Ths miracle of graphene is that non- 
relativistic electrons scattering against a lattice potential 
experience a low energy world in which, in non-trivial 
regards, they behave in ways reminiscent of Dirac's rel- 
ativistic fermions: the Klein paradox, Landau quantiza- 
tion, the fine structure constant, and so forth. An am- 
bitious but natural question to ponder is whether this 
coarse grained graphene world might mimic some aspects 
of gravitational structure? In Cartan's generalization of 
Einstein's geometrical formulation of gravity, torsion and 
curvature can be put in one-to-one relation with the dis- 
locations and disclinations, the topological defects of the 
crystal lattice [l| . The analogy is incomplete in the sense 
that crystal spaces are non-diffeomorphic. The general 
covariance of space-time translates into the requirement 
that the action of the medium should be independent of 
arbitrary elastic deformations, something obviously vi- 
olated by crystal and defected crystal geometry, which 
therefore correspond to spaces with a preferred metric or 
a 'fixed frame'. In addition, in the crystal the transla- 
tions and rotations of Galilean space are broken to dis- 
crete subgroups and this implies that both curvature and 
torsion arc quantized in units of the discrete Burgers and 
Frank vectors, the topological invariants of the disloca- 
tions and disclinations respectively. One can nonetheless 
still study the transport of matter in such topologically 
non-trivial fixed frame backgrounds; several such analo- 
gous gravity systems have been identified, including the 
sound waves of superfiuid He^ and the nodal Fermions 
of the He^ A phase, which percieve the hydrodynamical 
flow fields as geometrical (Christoffel) connection Q. 

To what extent does this analogy extend to the Dirac- 
like fermions ingraphene? We are inspired by the pre- 
vious work 0] demonstrating that the holonomy ac- 
cumulated by electrons in graphene encircling a disclina- 
tion (cone) coincides with that associated with a Dirac 
fermion encircling the conical singularity, the entity en- 
capsulating curvature in 2-1-1 dimensional gravity. How- 
ever, in order to complete the identification these earlier 
works added an ad-hoc U{1) gauge flux to the conical 
singularity, acting with opposite sign on the valley quan- 



tum numbers of the graphene electrons, raising the issue 
of whether the identification is merely coincidental. Here 
we will settle these matters by focussing on the influence 
of dislocations, corresponding to torsion in the gravita- 
tional analogy. Torsion is a less familiar aspect of the geo- 
metrical formulation of gravity [EH, [1,0]. Ignored at first 
by Einstein, it was introduced by Cartan [^ as an a priori 
ingredient of a geometrical theory. It was later pointed 
out by Kibble [vf that its inclusion becomes necessary 
in the presence of spinning particles, as their spin cur- 
rents source torsion in a dynamical space-time, though 
whether torsion propagates in the space-time is a matter 
to be settled by observation, an as yet open question be- 
cause torsional effects turn out to be too weak to be mea- 
sured with present day experimental technology, but the 
situation is different in the 'analogous' graphene system. 
Dislocations correspond with large 'fixed frame' localized 
torsion sources. We demonstrate that the holonomies as- 
sociated with graphene electrons encircling dislocations 
resemble those coming from the most natural implemen- 
tation of torsion in the connection of doubled fundamen- 
tal fermions, in the case that their cones would be dis- 
placed away from zero momentum. It is just the fact 
that the discreteness of 'graphene geometry' is remem- 
bered exclusively by the long wavelength fermion modes 
by the large momenta where the Dirac cones reside, and 
this is surely different from the way that Planck scale dis- 
creteness (when it exists) affects fundamental fermions. 
We subsequently show that the mysterious U{\) fiux of 
the graphene disclination has precisely the same origin, 
bringing us to the conclusion that the parallel transport 
of electrons in graphene with dislocations and disclina- 
tions is in the long wavelength limit identical to that 
of Dirac fermions living at large momenta in a 2-1- 1-d 
Cartan-Einstein spacetime with torsion and curvature. 

Torsion in Elasticity and its Coupling to Fermions - 
Within a geometric formulation of elasticity theory, dis- 
locations become sources of torsion (see [l|, l^, [iP^ and ref- 
erences therein) , stemming from their translational char- 
acter. Torsion T assigns a vector to an infinitesimal area 
element at each point in space, measuring the non-closure 



2 



of a loop obtained by parallel transport of the two in- 
finitesimals forming the "edges" of the given surface ele- 
ment along each other [l|, 13 • The definition makes 
this vector completely analogous to the Burgers vector in 
a crystal lattice. This gravity/geometry analogy has been 
verified in familiar electron systems, produ cing results 
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compatible with the tight-binding approach 
The metric is connected to the displacement field in the 
crystal via (we use 77^1/ to designate the flat Minkowski 
metric in a cartesian basis 77 — diag(l, — 1, — 1), and 
ds^ = g ^^dx^ dx^ for the invariant distance): 



gij = rUj + diUj + dji 



(1) 



Time is essentially decoupled from space in this con- 
densed matter system (g^o = 0). 

Formally, torsion is defined as a vector valued 2-form 
on space-time, = T/^^da;^ A dx" , with ^, e {0, 1, 2}, 
and A the wedge product of differential forms. In this 
work the relevant example is of one dislocation at the ori- 
gin of two dimensional space, with Burgers vector b, with 
a corresponding torsion (see [H, and references therein) 



6° 5{x) dx Ady . 



(2) 



The flux of this form through any area containing the ori- 
gin is given by the Burgers vector // T° = 6". As usual, 
the metric of Eq. ([T]) determines the symmetric part (in 
the lower indices) of the connection, i.e. the Christoffel 
piece, while torsion adds additional information about 
parallel transport in space, being related to the asym- 
metric part of the connection T^^ — ^(r^^ ~r^^). 
The geometry is consistently defined only if the Einstein- 
Cartan structure equations are satisfied: 



R = dT + -[T,T] 
T = d(3 + T AP , 



(3a) 
(3b) 



r^j^dx'^ is written as a matrix 



where the connection T 
valued 1-form, the curvature ii is a matrix valued 2-form, 
and p is a frame (i.e. p°'{x) are the dual basis vectors of 
the tangent space at x). 

Let us now define the 2-1-1 dimensional structure of the 
graphene Dirac equation 14 , 1^ by identifying the Dirac 



matrices as 7 — T3 (8) 173, 7 = irs (8) C72, and 7 — —its 
(Ti, which satisfy the Dirac algebra {7^,7^} = 2ri'^^. The 
r Pauli matrices act in the space of the valley index K± , 
while the a Pauli matrices act on the sublattice {A/B) 
degrees of freedom. 

Since spin is defined with respect to rotations acting 
in a tangent frame, to study the equation of motion of a 
spinning particle we must introduce Q an orthonormal 
set of basis vectors E"", connected to the holonomic frame 
dx^ through the vielbein (here dreibein) dx^ = e^i?° 
(and the inverse (ej^)~^ = e°). Intuitively since spinors 
are square roots of vectors, we need the square root of 



the metric to define spinors on curved manifolds, and the 
vielbeins provide this square root eer] ^ g; more formally 
they intertwine the representations of the orthogonal ro- 
tation group and its covering spin group. Then the rel- 
evant (zero mass) Dirac equation in a curved torsionful 
background takes the form 



with the covariant derivative given via 



9^ - -CJ^ab7"7'' 



(4) 



(5) 



The displacement field corresponding to a dislocation 
situated at the origin in two spatial dimensions is well 



known and via eq. ([T]) determines the metric 17 1 
(note that b is to be regarded as infinitesimal in the con- 
tinuum theory, so that we retain only linear terms). For 
simplicity we take the Poisson ratio cr = 0, and fix b to 
point along the x axis. The strategy is then to find an 
orthonormal basis E'^ on this space, i.e. 



(6) 



and then the spin connection from eq. (j3bp by using the 
physical input about the defect in ([2]). For the basis we 

get 



E^ 
E^ 



1 — sin 6 cos (j) 

7^- cos (b 1 — TT— sin f 

27rr ^ 27rr ' 



dr 



(7) 



and uj\ 



= ~a(p— {cos (f)dr + sin (f)d(f)). It is note- 
worthy that the matrix of 1-forms is antisymmetric as 
it represents the rotation of the orthogonal basis during 
parallel transport. The first term d4> appears simply due 
to the use of pol ar coordinates, and is responsible for a 



term — 7^/2r 15|, Il8| in Q, as well as for an unphysical 
delta function in the curvature obtained through (|3a|) . 
a byproduct of the singularity in the polar coordinate 
system. 

The spin connection produces a trivial holonomy for 
the Dirac spinor, but a non-trivial topological action is 
present in the basis (there is in fact some freedom 
in the Einstein-Cartan formalism to move torsion effects 
between the basis one-forms and the spin-connection, ob- 
vious in eq. pb)) ). The connection encodes for the inte- 
grable elastic deformation around the dislocation, and 
the vector field corresponding to the u) "potential" in ([5|) 
follows the deformation of the crystal due to the miss- 
ing row of atoms. More quantitavely, the singularity in 
displacement encoding the topological defect is fully con- 
tained in Up = —^\n{x + iy), due to ^ uo = —bx [3]. 
In this case a; = 0. Another exaniple is the elasti- 
cally unrelaxed displacement field in l*!, corresponding 
to E^ =dx + ^d(j), E'^ = dy, and = 0. 

The above examples are instructive in emphasizing the 
relation JJ T'' = JJ dE" = § E" = b°-, where one should 
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note that d{d4>) = 2TT6{r)S{(j))dr Ad(j) due to (j) not being a 
holonomic coordinate at the origin. Obviously, the topo- 
logically non-trivial part of E will always be in the form 
of a -^dcj) correction to the basis vector E"" along the 
Burgers direction, b ■ x"" — b. This can easily be checked 
explicitly for our setup in eq. ([7]), if the i?" basis is ro- 
tated to {E'^jEy}, with our previous choice of Burgers 
vector. The topological action of the dislocation on the 
Dirac electron should then be viewed as a Berry phase 
arising from the term 

= z7^9m exp (^J dx^fl^d^j = , (8) 

where = f^^i—li) is the perturbation proportional 

to the Burgers vector. The non-trivial holonomy (Berry 
phase) is responsible for the salient feature of long range 
influence of the crystal defect [3, [l^, [13] , taking the value 



H{b) 



(9) 



where we recognize the Volterra operation of translating 
the wavefunction by the Burgers vector to describe the 
topology of a dislocation. However, the correct holon- 
omy follows from the effect of translation by b (which is 
of order of a lattice constant) on the true Bloch wave- 



function [15|, |18| , in other words 



Hlattice 



ibKT3 



(10) 



The connection is striking and pleasing, because the con- 
tinuum translation generator — ?V is replaced by a trans- 
lation generator Kt^ of the underlying lattice wavefunc- 
tion, which is a finite momentum {K±) state. 

Eq. @ encapsulates the essence of arguments relating 
the vielbein and the gauge field of Poincare (here Eu- 
clidean) group translations in gauge theories of gravity, 
and one might consider a continuum-limit theory living 
on a background where the Euclidian group generators 
are modified to accommodate the lattice constant sized 
interactions of the defects and the finite momentum ref- 
erence (Fermi) state. 

Curvature and Disclinations ~ In the case of disclina- 
tions, the associated curvature exists in 2-1- 1-d as conical 
singularities, and has been considered in the graphene 
lattice [1, 2m, 13 ■ However, special care has to be taken 
to include the exchange of Fermi points, i.e. the internal 
degree of freedom, that occurs for specific opening angles, 
by using an additional gauge field with only r operator 
structure. Therefore an additional gauge field is intro- 
duced, alongside the curvature. Following the discussion 
in the previous section it becomes clear that it is more 
consistent to view the additional Fermi point effect as 
a change in the generator of rotations for the graphene 
Dirac spinor. 

The correct holonomies in the presence of a discli- 
nation with the fundamental opening angles at the 



origin, obtained by the Volterra construction, are [1, 
H, [111 H{2tt/3) = exp(-i^f-) and H{Tr/3) = 



-ZTi exp (— z^^). Note that rotating by tt/S maps the 
Fermi points into each other, hence the Ti matrix. We 
rewrite this in an illuminating way (0 is the angle of discli- 
nation): 



H{9 



,-'ie(o-3+3Ti)/2 



(11) 



where we see the spinor rotation (half-angle) generator 
CT3/2 replaced by ((73 -I- 3ri)/2, in order to accommo- 
date the finite lattice constant effect due to the exis- 
tence of two electron species, at finite momenta K±. 
This is a generalization to the spinor case of the obser- 
vation that the disclination holonomy is the represen- 
tation of the rotation operator by the defect opening 
angle [1^. It stems from the fact that the spin con- 
nection term, which produces the non-trivial holonomy 
in this case, is actually given by the rotation generator 



.[7^,7" 



OJ. 



Ml2 — 



£3. 
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and fixes the curva- 



ture 2-form i?2 = — i?f = duj = 96{x)dx A dy. Note that 
the uj^ matrix is again antisymmetric. 

General Torsion Couplings - Here we attempt to iden- 
tify additional possible couplings of torsion to the specific 
electronic degrees of freedom in graphene, based on gen- 
eral considerations (see [HI for a similar analysis in a 
different condensed matter system). 

The Riemann-Cartan curved space with torsion is de- 
fined by eq. ([3]), and fixed through the choice of the 
connection (once a tangent basis is specified), which it- 
self provides the covariant derivative to be used in the 
Dirac equation, eq. (j4|). This coupling of geometry to 
the spinor can in principle be extended by additional 
scalar terms containing torsion, which might follow from 
the choice of an action for the full gravity+matter the- 
ory [3, [l^, 17, 24 1, or in some cases only by an ad- hoc 
choice. These terms are linear in torsion at the least, 
and so effectively behave as a delta function potential in 
space (eq. ([2])). Obviously this makes no contribution 
to a holonomy, but the discussion is interesting from a 
general viewpoint. 

If we choose to start from a covariantized Dirac La- 
grangian in 2-1-1 dimensions (see [6| for the treatment 
of 3-1- 1-d), we get an additional term in the Dirac equa- 
tion 27° (Va + T^b) — (written in anholonomic co- 
ordinates, with the covariant derivative becoming Vq = 
da — jUJabcl^^'^)- At this poiut we can extract all similar 
torsion content from the covariant derivative , by sep- 
arating the antisymmetric part of the connection. Again 
in 2-t-l-d we get Va = Va - ^Tj - iTabdV, where Va 
contains only the Christoffel symbol part of the connec- 
tion. The Dirac equation reads 



1, 



1 



b-bcd. 



T, 



bed 



* = 0, 



(12) 



with the formally defined "traditional" 7^^ = 17^7^7^ 



4 



rp b 



Ts (8) 1. It seems that since the topological effect of the 
dislocation is present strictly in the E'^ basis, which stems 
from the singular displacement field through the metric 
(eqs. (dl and ©), it is enough to retain the Christoffel 
connection part of Va, as if there was no torsion (the 
additional terms in eq. do not contribute) . One must 
note, however, that torsion cannot be simply disregarded, 
as it is present in the space due to eq. (l3b|) . 

Our form of T (eq. ((21)) constrains the polar vector 
{zxb)aS{x) to be orthogonal to the Burgers vector, 
and this is the only possible polar term. Considering ax- 
ial vector couplings generally, in the relevant 2+1 dimen- 
sional case, there is no traditional 7^ matrix which is in- 
dependent of the 7° algebra, and which could be used to 
reduce the spinors to Weyl components, because it com- 
mutes, instead of anticommutes, with the 7°. However, 
in the case of graphene we are dealing with a reducible 
representation of the Clifford algebra, built out of two 
irreducible ones (one at each K± Fermi point). For this 
case, there exists a Jnew matrix, which can be defined for 
the present odd dimensional situation and having all the 
properties of acting in even dimensions ^25^1 . The Jnew 
represents the parity transformation which mixes the two 
irreducible representations, i.e. in our case it must map 
between and K spinor components (note that they 
are connected through parity, as — while in 

contrast the dislocation gauge coupling, which it should 
reproduce, acts via phase shifts without coupling the two 
K points, i.e. it is of the form. The above observa- 
tions do not prevent the appearance of terms containing 
7^ , and the last term in eq. (fT^ is of such a form, but 
it happens to be identically zero due to the contraction 

S'^'^Tabc = 0. 

To further connect with the lattice dislocation coupling 
eq. (jlOp . one could consider the generalization of form- 
ing scalars making use also of the K vector. The al- 
lowed combinations are e'''"'Tf,^^Kdj! and e'^'^'T^^Kcj^ , 
but neither is usable. The first one has the free index 
timelike a = (contributing a time dependent Berry 
phase constant in space), due to non-zero Tf^f having 
purely spacelike indices. The second term has the same 
feature {K also has no time component), although it has 
the correct matrix form e'^'"'Ti^fKc'yf = —b ■ Kt36{x). 

Conclusions - We have shown how electrons in defected 
graphene can be viewed as moving in a geometry with 
curvature and torsion, with all the topological lattice ef- 
fects included in an appropriate adjustment of the under- 
lying space symmetry generators. This is a fresh view on 
the subject in graphene, treating both types of defects 
equally, while matching them clearly with their govern- 
ing symmetry sectors. We hope this perspective will aid 
in understanding systems with many defects, where the 
fact that the holonomies are non-Ableian, renders them 
highly non-trivial. Such a treatment must necessarily 
deal with intricacies of the continuum limit of the lattice 
description, as we have emphasized. 
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